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Simulations of galactic disks including a dark baryonic 
component 



Y. Revaz^'^, D. Pfenniger^, F. Combes^, F. Bournaud^ 
Abstract. 

ACDM numerical simulations predict that the "missing baryons" reside 
in a Warm-Hot gas phase in the over-dense cosmic filaments. However, there 
are now several theoretical and observational arguments that support the fact 
that galactic disks may be more massive than usually thought, containing a 
substantial fraction of the "missing baryons" . Hereafter, we present new N- 
body simulations of galactic disks, where the gas content has been multiplied 
by a factor 5. The stability of the disk is ensured by assuming that the ISM 
is composed out of two partially coupled phases, a warm phase, corresponding 
the observed CO and HI gas and a cold coUisionless phase corresponding to the 
unseen baryons. 



1. Introduction 



While the ACDM scenario predicts that the dark matter in galactic disks is dis- 
tributed in a dynamically hot heavy halo, there are now several observational 
facts showing that a fraction of this rnatter may be b aryonic, lying in the o uter 
disk of galaxies. Long ago, iBosmal ( 19781 ) (see also iHoekstra et al.l [2OO1I ) ob- 
served a correlation between the HI and dark matter in spiral galaxies, indicating 
a physical link between these two components. This HI-DM c orrelation is con- 
firmed by the baryonic Tully-Fischer relation (jPfenniger &: Re vaz 2005). Dark 
matter in the disk of galaxies is supported by dynamical arguments based on the 
asymmetries of galaxies, either in the plane of the disks or transversal to it. Large 
spiral structure present in the very exten ded HI disk o f NG C 2915 is na turally 
expla ined by a quasi self-gravitating disk ( Bureau et aliri999l : iMasset h Bureau 
I2OO3I ) . Self-gravitating disks are also ver tically unst able and give a natu ral expla- 



nation to the old problematic of warps (iRevaz &: P fenniger 2004). Thi lker et al 



(|2005l 1 observed UV-bright sources in the extreme outer disk of M83, revealing 
unexpected star formation, indicating that molecular gas, must be present in 
abundance in the outer disk of galaxies. Other arguments supporting the idea 
that dark matte r could reside in the g alactic disk in form of cold molecular gas 
are discussed in lPfenniger et all ( 199 
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2. A new multiphase ISM model 

Adding dark baryons in a galactic disk poses the problem of the global sta- 
bility of the system. This has been solved by assuming that the ISM is com- 
posed out of two dynamically decoupled phases. A visible gas phase (vg) corre- 
sponding to the CO and HI observed gas, having velocity dispersions of about 
lOkm/ s, and a dark gas phase ( dg) of very cold, clumpy, quasi-collisionless 
phase (jPfen niger Combed 19941 ). This latter is dynamically warmer (higher 



velocity dispersions) ensuring the disk stability. The transition between the 
two phases depends on the gas temperature and could be computed using the 
first thermodynamic principle combined with ISM heating and cooling processes. 
However, this method strongly depends on the gas density that is only a poor 
spatial averaged quantity in numerical simulations. In order to avoid this den- 
sity dependence, we assume that the transition times between the two phases 
statistically depends simply on the local UV flux which is assumed to be the 
dominant ISM heating process. The following relations for the transition times 
has been chosen as follow : 

where x represents the normalized local UV flux. The parameter r gives the 
time scale of the transition, independently of the outer flux. At equilibrium, 
these latter equations provide simple relations between visible gas and dark gas 
density. Near young stars characterized by strong UV flux, the transition time 
between visible to dark gas is long and thus the ratio visible/dark gas is high. 
On the contrary, in the outer galactic regions dominated by a weak extragalactic 
UV flux, the transition between visible to dark gas is very short and the dark 
gas dominates there. 

3. Application to isolated disks 

In this section, we compare the evolution of two models of isolated spiral galaxies 
different by their dark matter content. The first model is a "classical" galaxy 
having a collisionless massive dark halo and no additional dark baryons. In 
the second model, a fraction of the dark halo mass is transferred into a dark 
baryonic gaseous disk (the baryonic gas is multiplied by a factor of 5) that takes 
part to the ISM cycle described earlier. The galaxy is constructed in order to 
share the same observational properties than the "classical" model and being 
consistent with Eq. [TJ a similar stellar component (exponential disk + bulge), 
a similar rotation curve (flat up to 40 kpc) , and a similar surface density of the 
visible gas (up to 20 kpc). The dark disk is as large as the dark matter halo, 
extending up to 90 kpc. Its contribution to the total surface density is negligible 
at the center, while it is large in the outer part, where the UV field is weak. The 
initial conditions have been computed by resolving the second moments of the 
Jeans equation. The radial velocity dispersion of the visible gas is set to lOkm/s 
while for the dark gas, it is set such that the Safronov-Toomre parameter Q to 
1, ensuring disk sta bility. The mod els have been run using a modified version of 
the Gadget-2 code (jSpringe]||2005l ). including the ISM cychng. Model 1 and 2 
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Figure 1. Comparison between the surface density of the visible gas of 
model 1 and visible and dark gas of model 2, after 2Gyr of evolution. The 
size of each box is 200 x 200 kpc. 



have respectively 524'285 and 804'751 particles ensuring a spatial resolution of 
250 pc. 

The evolution of the disks shows a similar global behavior between the two 
models. After 4Gyr both models have a very similar rotation curve and both 
have formed a comparable bar. Despite having multiplied the baryonic gaseous 
disk by a factor 5, model 2 is globally stable. It's stability is ensured by the 
larger velocity dispersion of the collisionless dark gas, being about twice the one 
of the visible gas. Fig. [T] compares the surface density of the visible disk after 
2 Gyr. While the azimuthal averaged surface density is similar between the two 
models, the spiral structure in the outer 20 kpc is very different. Model 1 has 
nearly no spiral structure while model 2 displays open and contrasted spiral 
arms up to the edge of the visible disk. The presence of the dark gas in model 
2 increases the surface density and makes the Safronov-Toomre parameter^ be 
constant around a value of 2. The self-gravity of the disk is increased making it 
quasi-stable. On the contrary, in model 1, the surface density is much smaller 
and Q is about twice larger, ensuring the stability of the disk. This behavior has 
to be related with puzzling observations of spiral struct ures in HI disk, diff icult 
to be explained by the self-gravity of the HI disk itself ( Bureau et al. 19991 ). 



4. Application to the missing mass in tidal dwarf problem 

Observation of the system NGC5291 revealed that the tidal dwarf galaxies 
present in the debris of NGC5291 have an unexpected dark matter content. 
In this section, we show that the presence of dark gas in the debris of galaxies is 
a natural product of the assumption that galaxies have an additi onal dark gas 
component. Using the initial condition of Bournaud et al.l ( 20071 ). we have re- 



produced the HI ring of NGC5291 using our new multiphase model. The initial 
spiral galaxy is pretty similar to the model 2 (see below), with a slightly more 
massive disk. The left and middle panel of Fig. [2] show the surface density of the 



^Here, the Safronov-Toomre parameter is computed by taking into account only the disk com- 
ponent : stars, visible and dark gas, for the surface density and radial velocity dispersion 
estimation. 



4 



gas, 400 Myr after the collision in a box of 200 x 200 kpc. Tidal dwarf galaxies 
have formed in a ring, out of the debris of the dissipative visible gas. The low 
velocity dispersion of the dark gas (compared to the dynamically hot dark halo) 
lets it fall into potential wells of the tidal dwarfs, adding a dark component to 
this latter objects. The typical rotation curve of the tidal dwarfs formed is dis- 
play on the right panel of Fig. [5J In these objects, the dark gas (dotted curve) 
plays the same role than a halo in a spiral galaxy. It allows the total rotation 
curve (plain curve) to stay straight up to 6 kpc, while the contribution of the 
visible baryons (star + gas) decreases. Nearly no dark halo mass is contained 
in these objects and its contribution to the total rotation curve (dashed-dotted 
curve) is negligible. In this latter figure, the upper error bars correspond to 
the observed velocity curve of NG C5291a and the latter t o the modeled veloc- 
ity curve of its visible component ( Bournaud et al.ll2007l ). A good agreement 



between model and observations is obtained. 
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Figure 2. Left: surface density of the visible gas, Middle: surface density of 
the dark gas. The size of each box is 200 x 200 kpc, Right: Velocity curve of 
a typical tidal dwarf galaxy, compared with data of iBournaud et al] (|2007f ). 



5. Conclusion 



New numerical simulations show that the hypothesis where galaxies are assumed 
to have an additional dark baryonic component is not in contradiction with 
galaxy properties. Multiplying the HI disk by a factor of 5 let the disk stable 
if the ISM is composed out of two dynamically decoupled phases. Such galaxy 
gives a natural solution to different puzzling problems like the the presence of 
spiral structures in the outer HI disks and the missing mass problem in collisional 
debris of galaxies. 
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